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chamber contour is made by considering the heat exchange between the com- 
bustion product flow in the thrust chamber and the coolant flow in the 
cooling jacket. The steady-state conditions that are considered require 
the temperatures of the combustion products, the chamber walls, and also 
the heat flux through the walls to remain constant at any point in time. 
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wall temperatures and the specific heat flux through the wall remain con- 
stant with time at any given point. 

1. The convective specific heat transfer rate on the hot gas-side 
is written 


where q = convective heat transfer rate per unit area, 

^wg 

T wg = hot 9 as ~ slde wal1 temperature, 

hg = hot gas heat transfer (film) coefficient, 

T = adiabatic wall temperature is that which would be 
wa attained by the surface of an adiabatic or insulating 
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axial distance from the throat = x, 

n _ = enhancement factor for associatedturbulence, surface 
E roughness of the tube, and curvature effects [13]. The 

accuracy of the enhancement factor significantly affects the 
heat transfer calculation and the resulting wall temperatures. 




The solution to the unknown q, T . T.„. T„ , and also h„ and h„ , can be 



a single wall and not of tubes, and ignoring the heat transfer through the 
walls which separate the coolant passage. 

The cooled surface area A. in contact with the hot gases of an arbi- 


trary section ax,, between the stations x. and x 



Figure 3. Up-pass coolant flow temperature analysis in a nozzle section 



The up-pass outlet coolant temperature of the ith section is 
calculated by 




ere the variables on the right hand side were described in equations 
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7, Convergence criteria for the iteration loop. 

A solution is obtained when the absolute value of the relative 
difference between the two total heat fluxes tEqn. 131 is less 
small value as shown: 
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computational sequence for a solution. 
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Case 1 : with Enhancement Factor, n E = 1.0 
Case 2 : with Enhancement Factor, nr as in Table 



oo 



cooling tube, total number of tubes, wall thickness, heat conductivity of 
the wall material, coolant flow rate, recovery factor, enhancement factor, 
initial coolant temperature distribution along the wall, initial gas-side 
wall temperature distribution, etc. are given in Table 1; whereas, Table 2 


lists the physical coolant properties as a function of temperature for an 
expected pressure range between 4500 psia and 6000 psia [13, 15] in order 
to determine the coolant flow heat transfer coefficient, the specific heat, 
thermal conductivity and viscosity at each x-station for a coolant tempera- 



scription of these input data are given in reference 12. 



The effects of curvature, associated turbulence, surface roughness. 
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ing specific heat transfer rate through the chamber wall is shown in Figure 
6. These results are also given in Table 5. The total heat flux through 
the wall for n£ = 1.0 is 4318.75 [Btu/Sec]. 

The computed temperature profiles for an enhancement factor distribu- 


tion as given in Table 1, are shown in Figure 7. The enhancement factor 
distribution was chosen such that the computed gas-side wall temperature 
profile matches the test data [16]. The specific heat flux for this case 
is shown in Figure 8. The total heat flux through the wall is 4633.96 



through the wall for this case, Q w = 4633.96 [Btu/Sec] in Figure 8 approxi- 
mates the projected test data well. 




chamber wall and the coolant flow heat absorption. This computational con- 
cept can be used with any boundary layer analysis program and is not re- 
stricted to a particular computer program to estimate regeneratively cooled 
thrust chamber temperature profiles mentioned above. As an application of 


this method, the boundary layer analysis computer program BLIMP was used 
to estimate the temperature profiles (Twg, Tw£, T£) of the regeneratively 
cooled RL10 rocket engine; the RL10 thrust chamber was chosen because of 
recently acquired test data. Since an emperical equation for the coolant 





showing the complexity of this type of analysis. 
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